Cryptic proteins translated from deletion-containing viral genomes dramatically expand the influenza virus proteome.

Productive infections by RNA viruses require faithful replication of the entire genome. Yet many RNA viruses also produce deletion-containing viral genomes (DelVGs), aberrant replication products with large internal deletions. DelVGs interfere with the replication of wild-type virus and their presence in patients is associated with better clinical outcomes. The DelVG RNA itself is hypothesized to confer this interfering activity. DelVGs antagonize replication by out-competing the full-length genome and triggering innate immune responses. Here, we identify an additionally inhibitory mechanism mediated by a new class of viral proteins encoded by DelVGs. We identified hundreds of cryptic viral proteins translated from DelVGs. These DelVG-encoded proteins (DPRs) include canonical viral proteins with large internal deletions, as well as proteins with novel C-termini translated from alternative reading frames. Many DPRs retain functional domains shared with their full-length counterparts, suggesting they may have activity during infection. Mechanistic studies of DPRs derived from the influenza virus protein PB2 showed that they poison replication of wild-type virus by acting as dominant-negative inhibitors of the viral polymerase. These findings reveal that DelVGs have a dual inhibitory mechanism, acting at both the RNA and protein level. They further show that DPRs have the potential to dramatically expand the functional proteomes of diverse RNA viruses.

Cryptic proteins translated from deletion-containing viral genomes dramatically expand the influenza virus proteome.

Productive infections by RNA viruses require faithful replication of the entire genome. Yet many RNA viruses also produce deletion-containing viral genomes (DelVGs), aberrant replication products with large internal deletions. DelVGs interfere with the replication of wild-type virus and their presence in patients is associated with better clinical outcomes as they. The DelVG RNA itself is hypothesized to confer this interfering activity. DelVGs antagonize replication by out-competing the full-length genome and triggering innate immune responses. Here, we identify an additionally inhibitory mechanism mediated by a new class of viral proteins encoded by DelVGs. We identified hundreds of cryptic viral proteins translated from DelVGs. These DelVG-encoded proteins (DPRs) include canonical viral proteins with large internal deletions, as well as proteins with novel C-termini translated from alternative reading frames. Many DPRs retain functional domains shared with their full-length counterparts, suggesting they may have activity during infection. Mechanistic studies of DPRs derived from the influenza virus protein PB2 showed that they poison replication of wild-type virus by acting as dominant-negative inhibitors of the viral polymerase. These findings reveal that DelVGs have a dual inhibitory mechanism, acting at both the RNA and protein level. They further show that DPRs have the potential to dramatically expand the functional proteomes of diverse RNA viruses.

Mutational signature dynamics indicate SARS-CoV-2's evolutionary capacity is driven by host antiviral molecules.

The COVID-19 pandemic has been characterised by sequential variant-specific waves shaped by viral, individual human and population factors. SARS-CoV-2 variants are defined by their unique combinations of mutations and there has been a clear adaptation to more efficient human infection since the emergence of this new human coronavirus in late 2019. Here, we use machine learning models to identify shared signatures, i.e., common underlying mutational processes and link these to the subset of mutations that define the variants of concern (VOCs). First, we examined the global SARS-CoV-2 genomes and associated metadata to determine how viral properties and public health measures have influenced the magnitude of waves, as measured by the number of infection cases, in different geographic locations using regression models. This analysis showed that, as expected, both public health measures and virus properties were associated with the waves of regional SARS-CoV-2 reported infection numbers and this impact varies geographically. We attribute this to intrinsic differences such as vaccine coverage, testing and sequencing capacity and the effectiveness of government stringency. To assess underlying evolutionary change, we used non-negative matrix factorisation and observed three distinct mutational signatures, unique in their substitution patterns and exposures from the SARS-CoV-2 genomes. Signatures 1, 2 and 3 were biased to C→T, T→C/A→G and G→T point mutations. We hypothesise assignments of these mutational signatures to the host antiviral molecules APOBEC, ADAR and ROS respectively. We observe a shift amidst the pandemic in relative mutational signature activity from predominantly Signature 1 changes to an increasingly high proportion of changes consistent with Signature 2. This could represent changes in how the virus and the host immune response interact and indicates how SARS-CoV-2 may continue to generate variation in the future. Linkage of the detected mutational signatures to the VOC-defining amino acids substitutions indicates the majority of SARS-CoV-2's evolutionary capacity is likely to be associated with the action of host antiviral molecules rather than virus replication errors.

SARS-CoV-2 in Domestic UK Cats from Alpha to Omicron: Swab Surveillance and Case Reports.

Although domestic cats are susceptible to infection with SARS-CoV-2, the role of the virus in causing feline disease is less well defined. We conducted a large-scale study to identify SARS-CoV-2 infections in UK pet cats, using active and passive surveillance. Remnant feline respiratory swab samples, submitted for other pathogen testing between May 2021 and February 2023, were screened using RT-qPCR. In addition, we appealed to veterinarians for swab samples from cats suspected of having clinical SARS-CoV-2 infections. Bespoke testing for SARS-CoV-2 neutralising antibodies was also performed, on request, in suspected cases. One RT-qPCR-positive cat was identified by active surveillance (1/549, 0.18%), during the Delta wave (1/175, 0.57%). Passive surveillance detected one cat infected with the Alpha variant, and two of ten cats tested RT-qPCR-positive during the Delta wave. No cats tested RT-qPCR-positive after the emergence of Omicron BA.1 and its descendants although 374 were tested by active and eleven by passive surveillance. We describe four cases of SARS-CoV-2 infection in pet cats, identified by RT-qPCR and/or serology, that presented with a range of clinical signs, as well as their SARS-CoV-2 genome sequences. These cases demonstrate that, although uncommon in cats, a variety of clinical signs can occur.

Direct Nanopore Sequencing of Human Cytomegalovirus Genomes from High-Viral-Load Clinical Samples.

Nanopore sequencing is becoming increasingly commonplace in clinical settings, particularly for diagnostic assessments and outbreak investigations, due to its portability, low cost, and ability to operate in near real-time. Although high sequencing error rates initially hampered the wider implementation of this technology, improvements have been made continually with each iteration of the sequencing hardware and base-calling software. Here, we assess the feasibility of using nanopore sequencing to determine the complete genomes of human cytomegalovirus (HCMV) in high-viral-load clinical samples without viral DNA enrichment, PCR amplification, or prior knowledge of the sequences. We utilised a hybrid bioinformatic approach that involved assembling the reads de novo, improving the consensus sequence by aligning reads to the best-matching genome from a collated set of published sequences, and polishing the improved consensus sequence. The final genomes from a urine sample and a lung sample, the former with an HCMV to human DNA load approximately 50 times greater than the latter, achieved 99.97 and 99.93% identity, respectively, to the benchmark genomes obtained independently by Illumina sequencing. Thus, we demonstrated that nanopore sequencing is capable of determining HCMV genomes directly from high-viral-load clinical samples with a high accuracy.

Adeno-associated virus 2 infection in children with non-A-E hepatitis.

An outbreak of acute hepatitis of unknown aetiology in children was reported in Scotland1 in April 2022 and has now been identified in 35 countries2. Several recent studies have suggested an association with human adenovirus with this outbreak, a virus not commonly associated with hepatitis. Here we report a detailed case-control investigation and find an association between adeno-associated virus 2 (AAV2) infection and host genetics in disease susceptibility. Using next-generation sequencing, PCR with reverse transcription, serology and in situ hybridization, we detected recent infection with AAV2 in plasma and liver samples in 26 out of 32 (81%) cases of hepatitis compared with 5 out of 74 (7%) of samples from unaffected individuals. Furthermore, AAV2 was detected within ballooned hepatocytes alongside a prominent T cell infiltrate in liver biopsy samples. In keeping with a CD4+ T-cell-mediated immune pathology, the human leukocyte antigen (HLA) class II HLA-DRB1*04:01 allele was identified in 25 out of 27 cases (93%) compared with a background frequency of 10 out of 64 (16%; P = 5.49 × 10-12). In summary, we report an outbreak of acute paediatric hepatitis associated with AAV2 infection (most likely acquired as a co-infection with human adenovirus that is usually required as a 'helper virus' to support AAV2 replication) and disease susceptibility related to HLA class II status.

Effects of culling vampire bats on the spatial spread and spillover of rabies virus.

Controlling pathogen circulation in wildlife reservoirs is notoriously challenging. In Latin America, vampire bats have been culled for decades in hopes of mitigating lethal rabies infections in humans and livestock. Whether culls reduce or exacerbate rabies transmission remains controversial. Using Bayesian state-space models, we show that a 2-year, spatially extensive bat cull in an area of exceptional rabies incidence in Peru failed to reduce spillover to livestock, despite reducing bat population density. Viral whole genome sequencing and phylogeographic analyses further demonstrated that culling before virus arrival slowed viral spatial spread, but reactive culling accelerated spread, suggesting that culling-induced changes in bat dispersal promoted viral invasions. Our findings question the core assumptions of density-dependent transmission and localized viral maintenance that underlie culling bats as a rabies prevention strategy and provide an epidemiological and evolutionary framework to understand the outcomes of interventions in complex wildlife disease systems.

Evaluating the potential of whole-genome sequencing for tracing transmission routes in experimental infections and natural outbreaks of bovine respiratory syncytial virus.

Bovine respiratory syncytial virus (BRSV) is a major cause of respiratory disease in cattle. Genomic sequencing can resolve phylogenetic relationships between virus populations, which can be used to infer transmission routes and potentially inform the design of biosecurity measures. Sequencing of short (<2000 nt) segments of the 15 000-nt BRSV genome has revealed geographic and temporal clustering of BRSV populations, but insufficient variation to distinguish viruses collected from herds infected close together in space and time. This study investigated the potential for whole-genome sequencing to reveal sufficient genomic variation for inferring transmission routes between herds. Next-generation sequencing (NGS) data were generated from experimental infections and from natural outbreaks in Jämtland and Uppsala counties in Sweden. Sufficient depth of coverage for analysis of consensus and sub-consensus sequence diversity was obtained from 47 to 20 samples respectively. Few (range: 0-6 polymorphisms across the six experiments) consensus-level polymorphisms were observed along experimental transmissions. A much higher level of diversity (146 polymorphic sites) was found among the consensus sequences from the outbreak samples. The majority (144/146) of polymorphisms were between rather than within counties, suggesting that consensus whole-genome sequences show insufficient spatial resolution for inferring direct transmission routes, but might allow identification of outbreak sources at the regional scale. By contrast, within-sample diversity was generally higher in the experimental than the outbreak samples. Analyses to infer known (experimental) and suspected (outbreak) transmission links from within-sample diversity data were uninformative. In conclusion, analysis of the whole-genome sequence of BRSV from experimental samples discriminated between circulating isolates from distant areas, but insufficient diversity was observed between closely related isolates to aid local transmission route inference.

The P681H Mutation in the Spike Glycoprotein of the Alpha Variant of SARS-CoV-2 Escapes IFITM Restriction and Is Necessary for Type I Interferon Resistance.

The appearance of new dominant variants of concern (VOC) of severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) threatens the global response to the coronavirus disease 2019 (COVID-19) pandemic. Of these, the alpha variant (also known as B.1.1.7), which appeared initially in the United Kingdom, became the dominant variant in much of Europe and North America in the first half of 2021. The spike (S) glycoprotein of alpha acquired seven mutations and two deletions compared to the ancestral virus, including the P681H mutation adjacent to the polybasic cleavage site, which has been suggested to enhance S cleavage. Here, we show that the alpha spike protein confers a level of resistance to beta interferon (IFN-ß) in human lung epithelial cells. This correlates with resistance to an entry restriction mediated by interferon-induced transmembrane protein 2 (IFITM2) and a pronounced infection enhancement by IFITM3. Furthermore, the P681H mutation is essential for resistance to IFN-ß and context-dependent resistance to IFITMs in the alpha S. P681H reduces dependence on endosomal cathepsins, consistent with enhanced cell surface entry. However, reversion of H681 does not reduce cleaved spike incorporation into particles, indicating that it exerts its effect on entry and IFN-ß downstream of furin cleavage. Overall, we suggest that, in addition to adaptive immune escape, mutations associated with VOC may well also confer a replication and/or transmission advantage through adaptation to resist innate immune mechanisms. IMPORTANCE Accumulating evidence suggests that variants of concern (VOC) of SARS-CoV-2 evolve to evade the human immune response, with much interest focused on mutations in the spike protein that escape from antibodies. However, resistance to the innate immune response is essential for efficient viral replication and transmission. Here, we show that the alpha (B.1.1.7) VOC of SARS-CoV-2 is substantially more resistant to type I interferons than the parental Wuhan-like virus. This correlates with resistance to the antiviral protein IFITM2 and enhancement by its paralogue IFITM3. The key determinant of this is a proline-to-histidine change at position 681 in S adjacent to the furin cleavage site, which in the context of the alpha spike modulates cell entry pathways of SARS-CoV-2. Reversion of the mutation is sufficient to restore interferon and IFITM2 sensitivity, highlighting the dynamic nature of the SARS CoV-2 as it adapts to both innate and adaptive immunity in the humans.

Evidence that genomic incompatibilities and other multilocus processes impact hybrid fitness in a rattlesnake hybrid zone.

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system.

Recent changes to virus taxonomy ratified by the International Committee on Taxonomy of Viruses (2022).

This article reports the changes to virus taxonomy approved and ratified by the International Committee on Taxonomy of Viruses (ICTV) in March 2022. The entire ICTV was invited to vote on 174 taxonomic proposals approved by the ICTV Executive Committee at its annual meeting in July 2021. All proposals were ratified by an absolute majority of the ICTV members. Of note, the Study Groups have started to implement the new rule for uniform virus species naming that became effective in 2021 and mandates the binomial 'Genus_name species_epithet' format with or without Latinization. As a result of this ratification, the names of 6,481 virus species (more than 60 percent of all species names currently recognized by ICTV) now follow this format.

SARS-CoV-2 Evolution and Patient Immunological History Shape the Breadth and Potency of Antibody-Mediated Immunity.

Since the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), humans have been exposed to distinct SARS-CoV-2 antigens, either by infection with different variants, and/or vaccination. Population immunity is thus highly heterogeneous, but the impact of such heterogeneity on the effectiveness and breadth of the antibody-mediated response is unclear. We measured antibody-mediated neutralization responses against SARS-CoV-2Wuhan, SARS-CoV-2α, SARS-CoV-2δ, and SARS-CoV-2ο pseudoviruses using sera from patients with distinct immunological histories, including naive, vaccinated, infected with SARS-CoV-2Wuhan, SARS-CoV-2α, or SARS-CoV-2δ, and vaccinated/infected individuals. We show that the breadth and potency of the antibody-mediated response is influenced by the number, the variant, and the nature (infection or vaccination) of exposures, and that individuals with mixed immunity acquired by vaccination and natural exposure exhibit the broadest and most potent responses. Our results suggest that the interplay between host immunity and SARS-CoV-2 evolution will shape the antigenicity and subsequent transmission dynamics of SARS-CoV-2, with important implications for future vaccine design.

Neutralizing antibodies provide protection against viruses and are generated because of vaccination or prior infections. The main target of anti-SARS-CoV-2 neutralizing antibodies is a protein called spike, which decorates the viral particle and mediates viral entry into cells. As SARS-CoV-2 evolves, mutations accumulate in the spike protein, allowing the virus to escape antibody-mediated immunity and decreasing vaccine effectiveness. Multiple SARS-CoV-2 variants have appeared since the start of the COVID-19 pandemic, causing various waves of infection through the population and infecting­in some cases­people that had been previously infected or vaccinated. Because the antibody response is highly specific, individuals infected with different variants are likely to have different repertoires of neutralizing antibodies. We studied the breadth and potency of the antibody-mediated response against different SARS-CoV-2 variants using sera from vaccinated people as well as from people infected with different variants. We show that potency of the antibody response against different SARS-CoV-2 variants depends on the particular variant that infected each person, the exposure type (infection or vaccination) and the number and order of exposures. Our study provides insight into the interplay between virus evolution and immunity, as well as important information for the development of better vaccination strategies.

Genomic epidemiology of SARS-CoV-2 in a university outbreak setting and implications for public health planning.

Whole genome sequencing of SARS-CoV-2 has occurred at an unprecedented scale, and can be exploited for characterising outbreak risks at the fine-scale needed to inform control strategies. One setting at continued risk of COVID-19 outbreaks are higher education institutions, associated with student movements at the start of term, close living conditions within residential halls, and high social contact rates. Here we analysed SARS-CoV-2 whole genome sequences in combination with epidemiological data to investigate a large cluster of student cases associated with University of Glasgow accommodation in autumn 2020, Scotland. We identified 519 student cases of SARS-CoV-2 infection associated with this large cluster through contact tracing data, with 30% sequencing coverage for further analysis. We estimated at least 11 independent introductions of SARS-CoV-2 into the student population, with four comprising the majority of detected cases and consistent with separate outbreaks. These four outbreaks were curtailed within a week following implementation of control measures. The impact of student infections on the local community was short-term despite an underlying increase in community infections. Our study highlights the need for context-specific information in the formation of public health policy for higher educational settings.

SARS-CoV-2 Omicron is an immune escape variant with an altered cell entry pathway.

Vaccines based on the spike protein of SARS-CoV-2 are a cornerstone of the public health response to COVID-19. The emergence of hypermutated, increasingly transmissible variants of concern (VOCs) threaten this strategy. Omicron (B.1.1.529), the fifth VOC to be described, harbours multiple amino acid mutations in spike, half of which lie within the receptor-binding domain. Here we demonstrate substantial evasion of neutralization by Omicron BA.1 and BA.2 variants in vitro using sera from individuals vaccinated with ChAdOx1, BNT162b2 and mRNA-1273. These data were mirrored by a substantial reduction in real-world vaccine effectiveness that was partially restored by booster vaccination. The Omicron variants BA.1 and BA.2 did not induce cell syncytia in vitro and favoured a TMPRSS2-independent endosomal entry pathway, these phenotypes mapping to distinct regions of the spike protein. Impaired cell fusion was determined by the receptor-binding domain, while endosomal entry mapped to the S2 domain. Such marked changes in antigenicity and replicative biology may underlie the rapid global spread and altered pathogenicity of the Omicron variant.

Tracking SARS-CoV-2 mutations and variants through the COG-UK-Mutation Explorer.

COG-UK Mutation Explorer (COG-UK-ME, https://sars2.cvr.gla.ac.uk/cog-uk/-last accessed date 16 March 2022) is a web resource that displays knowledge and analyses on SARS-CoV-2 virus genome mutations and variants circulating in the UK, with a focus on the observed amino acid replacements that have an antigenic role in the context of the human humoral and cellular immune response. This analysis is based on more than 2 million genome sequences (as of March 2022) for UK SARS-CoV-2 data held in the CLIMB-COVID centralised data environment. COG-UK-ME curates these data and displays analyses that are cross-referenced to experimental data collated from the primary literature. The aim is to track mutations of immunological importance that are accumulating in current variants of concern and variants of interest that could alter the neutralising activity of monoclonal antibodies (mAbs), convalescent sera, and vaccines. Changes in epitopes recognised by T cells, including those where reduced T cell binding has been demonstrated, are reported. Mutations that have been shown to confer SARS-CoV-2 resistance to antiviral drugs are also included. Using visualisation tools, COG-UK-ME also allows users to identify the emergence of variants carrying mutations that could decrease the neutralising activity of both mAbs present in therapeutic cocktails, e.g. Ronapreve. COG-UK-ME tracks changes in the frequency of combinations of mutations and brings together the curated literature on the impact of those mutations on various functional aspects of the virus and therapeutics. Given the unpredictable nature of SARS-CoV-2 as exemplified by yet another variant of concern, Omicron, continued surveillance of SARS-CoV-2 remains imperative to monitor virus evolution linked to the efficacy of therapeutics.

Selection Analysis Identifies Clusters of Unusual Mutational Changes in Omicron Lineage BA.1 That Likely Impact Spike Function.

Among the 30 nonsynonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (1) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (2) interactions of Spike with ACE2 receptors, and (3) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron overall previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

Quantifying and Cataloguing Unknown Sequences within Human Microbiomes.

Advances in genome sequencing technologies and lower costs have enabled the exploration of a multitude of known and novel environments and microbiomes. This has led to an exponential growth in the raw sequence data that are deposited in online repositories. Metagenomic and metatranscriptomic data sets are typically analysed with regard to a specific biological question. However, it is widely acknowledged that these data sets are comprised of a proportion of sequences that bear no similarity to any currently known biological sequence, and this so-called "dark matter" is often excluded from downstream analyses. In this study, a systematic framework was developed to assemble, identify, and measure the proportion of unknown sequences present in distinct human microbiomes. This framework was applied to 40 distinct studies, comprising 963 samples, and covering 10 different human microbiomes including fecal, oral, lung, skin, and circulatory system microbiomes. We found that while the human microbiome is one of the most extensively studied, on average 2% of assembled sequences have not yet been taxonomically defined. However, this proportion varied extensively among different microbiomes and was as high as 25% for skin and oral microbiomes that have more interactions with the environment. A rate of taxonomic characterization of 1.64% of unknown sequences being characterized per month was calculated from these taxonomically unknown sequences discovered in this study. A cross-study comparison led to the identification of similar unknown sequences in different samples and/or microbiomes. Both our computational framework and the novel unknown sequences produced are publicly available for future cross-referencing. Our approach led to the discovery of several novel viral genomes that bear no similarity to sequences in the public databases. Some of these are widespread as they have been found in different microbiomes and studies. Hence, our study illustrates how the systematic characterization of unknown sequences can help the discovery of novel microbes, and we call on the research community to systematically collate and share the unknown sequences from metagenomic studies to increase the rate at which the unknown sequence space can be classified.

Intra-host analysis of hepaciviral glycoprotein evolution reveals signatures associated with viral persistence and clearance.

Even 30 years after the discovery of the hepatitis C virus (HCV) in humans there is still no vaccine available. Reasons for this include the high mutation rate of HCV, which allows the virus to escape immune recognition and the absence of an immunocompetent animal model for vaccine development. Phylogenetically distinct hepaciviruses (genus Hepacivirus, family Flaviviridae) have been isolated from diverse species, each with a narrow host range: the equine hepacivirus (EqHV) is the closest known relative of HCV. In this study, we used amplicon-based deep-sequencing to investigate the viral intra-host population composition of the genomic regions encoding the surface glycoproteins E1 and E2. Patterns of E1E2 substitutional evolution were compared in longitudinally sampled EqHV-positive sera of naturally and experimentally infected horses and HCV-positive patients. Intra-host virus diversity was higher in chronically than in acutely infected horses, a pattern which was similar in the HCV-infected patients. However, overall glycoprotein variability was higher in HCV compared to EqHV. Additionally, selection pressure in HCV populations was higher, especially within the N-terminal region of E2, corresponding to the hypervariable region 1 (HVR1) in HCV. An alignment of glycoprotein sequences from diverse hepaciviruses identified the HVR1 as a unique characteristic of HCV: hepaciviruses from non-human species lack this region. Together, these data indicate that EqHV infection of horses could represent a powerful surrogate animal model to gain insights into hepaciviral evolution and HCVs HVR1-mediated immune evasion strategy.